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ABSTRACT 


Thermal, convection, and rotational fluid flow analyses were conducted 
of a prototype of the SPAR experiment 76-36. A constitutional supercooling 
criterion (CSC) was calculated from the thermal data. The study results 
include the following: 


Cooling rates have been identified for generating 
thermal conditions such that freezing occurs 
within the allocated low-g time; thermal and solutal 
convection velocities are of the order of lO”^ cm/sec 
at 10"^g; and the CSCs are approximately those of 
the SPAR experiment 74-21. 


• An experiment orientation with respect to the rocket 
spin axis has been identified which should ensure 
optimum damping times of rotational motion. 
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I. Objectives 

The objectives of this study are: 


A. Provide analytical support for fluid flow phenomena in SPAR 
Experiment 76-36. 

B. Provide comparison analyses between SPAR Experiments 
74-21 and 76-36 for the purpose of experiment design. 


II. Statement of W ork 


The contractor shall provide the necessary personnel and facilities to 
perform the following tasks: 


Task 1: The contractor shall map the thermal profile in the 

melt for each of the three alloy systems in SPAR 
Experiment 76-36. 

Task II: The contractor shall map anticipated fluid flow due 

to solutal and buoyancy driven convection in the 
sample materials during the Sounding Rocket Flight. 

Task III: The contractor shall compare the low -gravity fluid 
flow obtained in SPAR Experiment 74-21 with that 
expected in SPAR Experiment 76-36. Design criteria 
for sizing the alloy samples will be provided so that 
the two SPAR Experiments will have identical initial 
fluid flow and thermal fields. 
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INTRODUCTION 

One of the first studies of dendrite crystal growth in low-g environments 
was the series of experiinents performed in the period 1975-77 on Black Brandt 
VC sounding rocket flights (SPAR Experiment 74-21). In these early experi- 
ments, dendrite crystals of ammonium chloride were grown from aqueous solu- 
tions cooled below supersaturation. The low-g crystallizations were recorded 
on film and subsequently compared with crystals crystallized on the ground. 

The low-g crystallization showed a number of different behaviors from those 
observed in the 1-g cases (Ref. 1). 

Two aspects of the ammonium chloride crystallizations in low-g are of 
interest. The basic phenomenon of dendritic crystallization in the absence of 
convection is of intrinsic importance. Of particular interest is the origin of 
the secondary nucleation of detached crystals occuring after the primary den- 
drite stalks have been nucleated and grown awhile. Some theories hold that 
convection breaks off delicate dendrite arms and sweeps them into the bulk 
of the melt. Other theories hold that the nucleation arises from an interplay 
of thermal and constitutional conditions. In a more applied vein, the crystal- 
lization of ammonium chloride can serve as a useful model for metal alloy 
crystallization both in 1-g and low-g. In order to be a useful model for metal 
alloy crystallization, however, convective fluid flow (velocities and flow patterns) 
and constitutional supercooling conditions need to be made similar. 

As discussed in a prior report (Ref. 2), making similar both convective 

and constitutional supercooling conditions for the case of 1-g gravity would 

involve a more lenghty study than can be undertaken at this time. In low-g 

conditions, however, convective velocities can be reduced to insignificant 
-4 -5 

levels (-^10 or 10 cm/sec). The only concern thus becomes that of 
similarizing constitutional supercooling conditions. Accordingly one of the 
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objectives of a second series of rocket experiments of dendrite growth is to 
freeze molten metal systems in low-g under themal conditions such that the 
constitutional supercooling conditions are generally the same as for the first 
ammonium chloride low>g crystallizations (SPAR Experiment 74-21). The 
data from the 74-21 and 76-36 should then show that ammonium chloride 
crystallizations can be a useful model for metal alloy crystallization as well 
as providing basic data on dendrite crystallization in general. 


The purpose of the present study is to provide prelimine'>*y thermal and 
convective analysis data for the 76-36 experiment in order to ensure that the 
experimental furnace will provide the desired constitutional supercooling and 
convective conditions. The approach used consists of four parts: 


• Thermal analysis of a preliminary furnace design to ensure 
that the metal samples are completely frozen in the avail- 
able low-g time period (about 5 to 7 min) of the SPAR rocket 
test, to provide data for determining the supercooling param- 
eter, and to provide thermal boundary conditions for the con- 
vection analyses. 

• Calculations of a supercooling criterion (C/RA ratio, where 
G is the temperature gradient at the crystallizing interface, 

R is the crystal growth rate, and A is a constant for a given 
solute concentration), for the various cases. 

• A convection analysis of the various cases to ensure that 
convective velocities will not exceed about 10"^ cm/sec in 
the low-g tests. 

• An analysis of damping times for fluid flow generated by 
rocket spin-up and spin-down. 


Each of these areas is discussed in detail in following sections as is a 
description of the general nature of the apparatus, materials, and procedures. 
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EXPERIMENT DESCRIPTION, APPARATUS, AND MATERIALS 


The SPAR 76-36 Experiment entitled **Comparative Alloy Solidification*' 
will consist of melting a metal alloy sample on the ground and bringing it up 
to some designated constant temperature (soak temperature) just prior to the 
rocket flight. During launch and about 2 min after the rocket enters the low-g 
period, the soak temperature will be maintained. Cooling at the sides of the 
cell will commence after the 2 min damping period* Solidification during the 
test will be somewhat as shown in the following two-dimensional representation* 


Cooling 


iUtaUtUdilia/iiL 


^—Adiabatic Top, Bottom, 
^ Front, and Back Walls 



Cooling 


Dendrites 


Fig. 1 - Two-Dimensional Schematic of Experiment 76-36 Crystallization 


The sample will be completely frozen by the time the low-g period is 
completed and the rocket starts its descent to earth. The frozen samples 
will be retrieved and analyzed on earth by various metallographic techniques. 

A prototype furnace design to accomplish the required soUdifications 
is shown in Fig. 2. 
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The crucible will be constructed of graphite and will have the dimensions 
(first-cut) shown in Fig. 3. 



X 

0.953 cm 


Graphite 



Top View of Crucible 


5.08 cm 


Liquid Cavity Dimensions 


Fig. 3 - Dimensions of Graphite Crucible 


The metal alloy systems to be used in the tests are given in Table 1. 
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Liquid and solid physical property data for the metal alloy systems as 
well as for the aqueous ammonium chloride solution of Experiment 74>21 are 
given in Tables 2 and 3. 

The following symbol identification is used in the tables. 

P = density 
c = heat capacity 
X s heat of fusion 
k' = thermal conductivity 
p = viscosity 
3^ = thermal expansivity 
3^ = solutal expansivity 
1 / = kinematic viscosity 

a = thermal diffusivity 
D = solutal diffusivity 
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THERMAL ANALYSES 

Temperature profiles as a function of time * ithin the cell at various 
conditions of cooldown were computer calculated. The Lockheed Thermal 
Analyzer Program was utilized. This program is applicable for cases of heat 
conduction within a body. Convection on the boundaries is handled by means 
of a heat transfer coefficient. The latent heat generated during solidification 
is handled in the program by ordering the program to decrease the heat capacity 
(sufficient to allow liberation cf the latent heat) at a node when the temperature 
reaches the freezing point. Such a procedure is fairly standard for computer 
calculations involving phase chauige. The cell was modeled according to the 
schematic shown in Fig. 4. 


Order of magnitude cooling rates were dete^-mined from a simple calcu- 
lation of cooling rate required to freeze each melt within 360 sec. The formula 
used for .his calculation is given by Eq. (1). 

Q = ^ (1) 


where 


thermal conductivity 

« T. 


initial 
360 sec 


freeze 


cooling rate 
thermal diffusivity 


The results of applying this equation to each of the melts are shown below. 



Unit 

Q (cal/^cm^-sec) 

AT (C) 

1 ^ 

Sn - 50 Pb 

0.50 

90 

i 

i 

j 

Al - 4.5 Cu 

2.86 

155 

i 

1 

[ , Sn - 3 Bi 

0.50 

113 
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Top View 



288 Node Points 
720 Conduction Paths 
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Plane of 
Symmetry 
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Node Points 
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Heater Power • 


Plane of Symmetry 


[Coolant 


Side View 

Furnace Canister 
(Full Size) 


Sample 


q (Cooling) 


Graphite 


^ Plane of 
Symmetry 


‘Node Points 
(288 Total) 


Fig. 4 * Computer Model of 76-36 Experiment 
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Using a value near the above estimate for Q, the Sn-50 Pb system tem- 
perature profiles are shown in Fig. 5. The model currently assumes that the 
ceramic heaters behave thermally like graphite and that all surfaces are 
adiabatic except those cooled by the gas. 

A number of cases were run varying soak temperature and cooling rates 
to provide data for calculation of the supercooling criterion, as discussed in 
the next section, and to provide boundary conditions for the convection analysis 
discussed in a following section. 

Isotherm profiles (Figs. 6 through 9) calculated by the Thermal Analyzer 
Program revealed a potentially significant furnace design factor. This factor 
consists of the sharply curved isotherms (isotherm profiles shown as views 
from the top looking down — "horizontal cuts"). The intent of the experiment 
is to freeze the metal with flat isotherms. These isotherms will be parallel 
to the cooling walls. Thus the final furnace design should take these curved 
isotherms into account, and design alternatives should be explored which will 
minimise these undesirable effects. 







Botherm Profiles £oi 
:er Quench Begins, ' 
t; Soak Temperature 








Fig. 7 - Furnace Isotherm Profiles for Sn Pb (Case 23, 

120 sec After Quench Begins, Values in Degrees 
Fahrenheit: Soak Temperature, 401 F (205 C)) 





Sample 


t M t M 1 



<0 u ^ 
C C sO 

feO - 



CONSTITUTIONAL SUPERCOOLING CRITERION (CSC) 


The liquidus temperature gradient (9T./8x' ) ahead of a solidifying 
interface can be shown to be given by the following equation (Ref. 3 gives a 
simple derivation). 


OTj^/ax' )^^_^ = - m R C° (1 - k)/kD 


where 

Tl = liquidus temperature 

= initial solute concentration 
m - slope of the liquidus curve 
k = segregation coefficient 
D = solute diffusivity 
R = rate of freezing 

X* = location of the moving solid/ liquid interface 

If the actual temperature gradient at the interface is G» we can specify 
a constitutional supercooling parameter (CSC) as follows: 

G G G 

(aT/8x' )^,^Q m R C° ( 1 - k)/ k D ^ 

The CSC parameter can thus be viewed as a dimensionless supercooling 
parameter. If CSC is less than one then the liquidus temperature gradient 
is steeper than the actual temperature gradient and constitutional super- 
cooling prevails. A value greater than one has the opposite meaning — the 
liquid melt is constitutionally superheated. 

Although it would seem that the CSC parameter — giving as it does the 
relative degree of supercooling — should also be an excellent indicator of 
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crystal morphology, no one has yet thought to test such a relationship* A 
number of prior studies have tried to correlate G/R or with crystal 

morphology (Refs. 4, 5 and 6) and have achieved varying success. One possible 
complication not considered thus far is that relatively slow crystal growth 
kinetics can result in considerable thermal as well as constitutional under- 
cooling. A test of the correlation of the CSC parameter with crystal morpho- 
logy is a study planned for the future by the present investigators. Tables 4 
through 6 show the value of the CSC parameter for various cases. As indicated 
in the tables, the CSC value at the time of initial freezing for the NH^CI in Ex- 
periment 74-21 was approximately 0.01. Thus, the results for all three metal 
systems of Experiment 76-36 are of the same order of magnitude as occurred 
in Experiment 74-21. Thus, the metal-model results of Experiment 74-21 
should be especially appropriate to those of Experiment 76-36. A more de- 
tailed account of the CSC throughout the whole period of freezing, however, is 
indicated in any future analysis. 
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CONVECTION ANALYSIS 

A convection analysis was conducted of two-dimensional, transient 

thermally-driven and solutally-driven convection for each of the three 

.4 

metal systems at a gravity level of 10 g . 



The container wall temperatures versus time, for each metal ca? 
studied, are shown in Figs. 10, 11 and 12. These temperatures were taken 
from the thermal analyses results. The Grashof and Raleigh numbers 
shown in these figures are "normalized," and the actual values of Gr and 
Ra for a given case are obtained by multiplying by AT/T^ (degrees Kelvin) 
or by AC/C^ (wt%) for thermal and solutal, respectively. 

The thermally driven convection results are shown in Figs. 13 through 
18. The results are in the form of velocity vector maps and isotherm plots 
at a given instant of time. For each metal, results are shown for an early 
time and a late time during the freezing cycle. The effects of the moving 
crystallization front on the flow were ignored, thus predicted flow velocities 


for later times will be higher than actual. Predicted trends, however, should 
remain valid. Maximum predicted velocities are of the order-of-magnitude 
10 ^ cm/ sec. 



The solutally driven convection results are shown in Figs. 19 through 
24. The data are presented in the form of velocity vector and isoconcentration 
plots for early and late times (in the freezing cycle) for each metal system. 
Again, maximum predicted velocities are less than 10*^ cm/sec. 














Isotherms 



Fig. 13 - Isotherms and Velocity Profiles at t=80 sec for SnPb Case 20 
(Thermal Convection): V = 4.7 cm/sec ; Isotherms in 

Degrees Celsius 







Isotherms 



Velocity Vectors 


' J 

' t t I / M 

^ t f M 

* t ! * 

♦ t \ ♦ 

‘ / f t N : j 


-V > ^ 


Fig. 14 - Isotherms and Velocity Profiles at t= 120 sec for SnPb Case 20 
(Thermal Convection); V = 5.2 ‘lO"^ cm/sec; Isotherms in 


Degrees Celsius 
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Fig. 16 > Isotherms and Velocity Profiles at t= 100 sec for Al Cu Case 12 

(Thermal Convection); Max. Velocity = 7.6 • 10"^ cm/sec; Isotherms 
in Degrees Celsius 
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Isoconcentration 


Velocity Vectors 



Fig. 24 - Isoconcentration and Velocity Profiles at t= 100 sec for SnBi 
Case 29 (Solutal Convection): V__„ = 1.9 *10"^ cm/sec; Iso- 


^ M / \i«rwav4 vca a 'w v c / f v rT1.3X 

concentrations in Wt. % Bismuth 







FLUID FLOWS GENERATED BY SPIN-UP AND DECAY TIMES 


During the launch phase of a SPAR rocket flight, various accelerations 
occur. These are summarized as follows: 


Stabilizing spin around longitudinal axis 240 rpm 

Lateral accelerations due to aerodynamic 4 5 g, 0-4 Hz 
forces (g-jitter) 

Maximum g due to rocket acceleration l6 g 


After about one minute after ignition, a "yo-yo" despin mechanism is deployed 

and the spin is brought down to about 0.2 deg/sec. Aerodynamic forces also 

.4 

become unimportant. Thus, the gravity environment becomes about 10 g at 
all points some 75 to 80 sec after ignition. 

In the 76-36 experiment molten metal at its soak temperature will be 
sent up and will experience two periods of rotational acceleration — a spin-up 
and a spin-down — in rapid order before entering the low-g period. The question 
of what fluid motion is generated by the spin-up and how long it takes to decay 
after spin-down thus becomes of crucial importance. It may be mentioned that 
thrust acceleration and g-jitter induced flow should be insignificant because 
the melts will be essentially isothermal, homogeneous fluids. 

The transient type of motion generated when a fluid filled cylinder is 
suddenly spun up is shown in Fig. 25. In Fig. 25 it should be kept in mind that 
the fluid sidewall layers are rotating while the center core remains irrotational 
for some time. / condition is finally reached where the entire body of fluid 
rotates as a unit. In the present simple analysis, it will be assumed that spin-up 
is complete in the liquid by the time the low-g phase of the flight is entered. 
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Fig. 25 - Schematic Drawing of Boundary Layers (Dashed) and Meridional 
Circulation Pattern Solid) for the Spin-Up of a Homogeneous Fluid 
in a Circular Cylinder (Ref. 9) 

A prior treatment of spin-down times (Ref. 7) indicates that appropriate 
spin-down times can be calculated according to one of two analytical treat- 
ments in the literature, i.e., 

McLeod's analysis (Ref. 8) if L/2R * 1 

Benton's analysis (Ref. 9) if L/2R « 1 

where E is the Ekman number (= v/il \J‘), v is the kinematic viscosity, il is 
the angular spin velocity, and L is the cylinder length, and R is the cylinder 
radius. 

The cell of the 76-36 experiment is a thin rectangle which could be 
placed in one of two positions with respect to spin axis of the rocket, i.e.. 

Fig. 26. 

1/4 

Values of E ' L/2R can thus be calculated for three cases, i.e., 

Table 7. 
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Fig. 26 - Two Possible Positions of 76-36 Cell with Respect to Rocket Spin Axis 


Table 7 

APPLICABLE THEORETICAL TREATMENTS FOR VARIOUS CASES 


L/2R 

Applicable 

Treatment 


6. 1 cm 

6.1 cm 

0.95 cm 

0.475 cm 

2.55 cm 

2.55 cm 

(Effective 

Radius) 



0.62 

0.057 

0.022 

McLeod 

Benton 

Benton 


The values given in Table 8 show the values of sp’a-down times for 

various cases and various treatments assuming a spin-down from about 240 

- 3 2 

to 1 rpm. A representative v value of 5 x 10 cm /sec was used in the 
calculations. 
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Table 8 

THEORETICAL SPIN -DOWN TIMES 


Case 


Theory 

Cyl. Radius, 
R (cm) 

Cyl. Length, 
L (cm) 

Time 

(sec) 

McLeod 

0.475 

Not Relevant 

25.0 

McLeod 

2.55 

Not Relevant 

80.5 

Benton 

2.55 

6.1 

52.5 

Benton 

2.55 

0.95 

8.1 


The values shown in Table 8 indicate that the shortest spin-down times 
occur in either Case (b), Fig. 26, or Case (a) if the assumption of an effective 
radius of 0.475 cm is valid, as it probably is. 


Further information on spin-down times was obU.ined from Dr. L. L. Lacy 
of NASA-MSFC. Dr. Lacy in the past performed a number of rotation experi- 
ments with different fluids (including mercury) and containers. He found he 
could correlate his results by the following simple formula 


T = |3 dyv 

where T is the damping time (for the velocity to reach 5 x 10 cm/sec), p 
is a constant and equal to 0.102 for an end velocity 5 x 10 ^ cm/sec, d is the 
cylinder diameter, and v is the kinematic viscosity. Damping times calculated 
by Dr. Lacy's formula are as follows: 


Diameter, 
d (cm) 

0.95 (Effective) 
5.10 


Time, 
T (sec) 

18.4 

530.6 


These results would indicate that Case (b). Fig. 26, would be u.ndesirable. 
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Dr. Lacy performed one experiment in which a rectangular cell filled 
with water was spun slightly off axis, i.e.. Fig. 27. It can be seen that the 
decay times are fairly short. Dr. Lacy also found that wetting and non-wetting 
conditions made a difference in the damping times (Table 9). 

The conflicting results of theory and experiment is not apparent at 
present, but the length-to-diameter ratio (aspect ratio) is undoubtedly the 
crucial factor. Dr. Lacy's results are probably more applicable to tall, fat 
containers than to short, thin containers. In light of Dr. Lacy's experimental 
results on off-center placement and available theory, an orientation such as 
shown in Fig. 28 appears advisable. Also, as long a damping period as pos- 
sible should be utilized, with 2 min as the absolute minimum. 
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Fig. 27 - Vortex Motion and Decay Time for Rectangular Container of Water 
for SPA Rocket Experiment (Thickness; 1.3 cm) (Courtesy of Dr. 
D. L. Dacy^ k^SFC) 
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Fig. 28 - Recommended Orientation of Experiment Cell with Respect 
to Rocket Spin Axis 


Table 9 


TOTAL DECAY TIME AS A FUNCTION OF VARIOUS P/kRAMETERS 
(Courtesy of Dr. L. L. Lacy, MSFC) 


Fluid 

Container 

Material 

Diameter 

and 

Length (cm) 

Wetting oi 
Nonwetting 

Decay 

Time 

(sec) 

Conditions 

Water 

Plastic 


W 

36 + 2 

Center Vial 


Plastic 

2.15 

5.6 

W 

31+3 

Vial 2" Off Center 
of Rotation 


Plastic 

2.15 

5.6 

W 

2 

Vial with Vertical 
Baffle 


Plastic 

1.15 

W 

10 +1 

Vial 2” Off Center 


Glass 

1.69 

5.0 

W 

20+1 

Center Vial 


Glass 

mmEmm 

N 

18.6 + 0.7 

Same as Above Except 
Coated with Krylon 


Glass 


N 

18.3 + 1.5 

Coated with Dow 
Corning 200 Fluid 


Glass 

1.69 

5.0 

N 

17.5 + 1.5 

Coated with RAM Mold 
Release 225 


Glass 

2.38 1 

7.4 



W 

46.8+0.9 

Center Vial 

Carbon 

Tetra- 

chloride 

Glass 


W 

30.3 + 0.9 

Center Vial 

Cyclo- 
hexanol 1 

Glass 

1.69 

W 

< 1 

Center Vial 

Mercury 

Plastic 

i 

! 

2.15 

2.8 

N 

186+10 

Center Vial 


Plastic 

1.69 

2.95 

N 

95 +5 

Center Vial 
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